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Abstract. We have explored the emission line trends in 
the integrated spectra of normal spiral galaxies of the 
Nearby Field Galaxy Survey (Jansen et al. 2000a and b), 
in order to investigate the relationships between dust ex- 
tinction, metallicity and some macroscopic properties of 
spiral galaxies. 

We found a very strong correlation between the H/3 
and Ha equivalent widths, implying that the difference 
between the extinction of the stellar and the nebular light 
depends only on the intrinsic colours of the galaxies, being 
larger for redder galaxies. 

The usual metallicity indicator for giant H II regions 
([O in] A4959,5007 + [O II] A3726,3729) / H/3 is not ap- 
propriate for integrated spectra of spiral galaxies, proba- 
bly due to metallicity gradients. Much better qualitative 
metallicity indicators are found to be [N n] A6584/[0 n] 
A3726,3729 and [N n] A6584/Ha, the latter having the 
advantage of being independent of reddening and being 
applicable also for galaxies with weak emission lines. 

With these indicators, we find that the nebular extinc- 
tion as derived from the Balmer decrement strongly cor- 
relates with the effective metallicity of the emission line 
regions. 

The overall metallicity of the emission line regions is 
much better correlated with galaxy colours than with mor- 
phological types. 

A Principal Component Analysis on a 7-D parame- 
ter space showed that the variance is produced, in first 
place, by the metallicity and parameters linked to the 
stellar populations, and, in second place, by the surface 
brightness, which is linked to the dynamical history of the 
galaxies. The absolute magnitude, related to the mass of 
the galaxy, comes only in the third place. 
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1. Introduction 

The interpretation of integrated spectra of nearby galax- 
ies becomes increasingly important with the profusion of 
studies of galaxies in clusters (e.g., Caldwell et al., 1993; 
Dressier et al. 1999; Poggianti et al. 1999; Balogh et al. 
1999) and galaxies at large redshifts (e.g., Colless et al. 
1990; Hammer et al. 1997; Lilly et al. 1998, Guzman et al. 
1998). 

Indeed, for the studies of distant galaxies, the spa- 
tial resolution obtained with ground based instruments 
makes any morphological classification difficult, and even 
from space, morphological studies are bound to be crude. 
Therefore, most information on distant galaxies will come 
from photometry or spectroscopy, and the sampled regions 
will cover a significant part - if not all - of the galaxy sur- 
face. 

Techniques have been devised to classify galaxies us- 
ing their spectra only (Morgan & Mayall; 1957, Sodre & 
Cuevas; 1994, 1997; Connolly et al. 1995; Zaritsky, Kenni- 
cutt & Huchra 1994; Folkes et al. 1996; Galaz & de Lap- 
parent 1998; Sodre & Stasihska 1999). It has been shown 
that the results correlate rather well with the galaxy 
morphological types. The spectral classification actually 
presents some advantages even for the study of nearby 
galaxies: whatever criteria are adopted, they allow an ob- 
jective and continuous classification of the galaxies, in 
contrast with the conventional morphological classifica- 
tion which is by nature subjective and provides discrete 
classes. In addition, any interpretation in terms of stellar 
content and evolution is directly related with the galaxy 
spectra, while the link with galaxy morphology is rather 
elusive. 

In a former paper (Sodre & Stasihska 1999, hereinafter 
SS99), based on the integrated spectra for normal spiral 
galaxies from the Kennicutt (1992) atlas, we have shown 
that emission line ratios and equivalent widths are im- 
pressively well correlated with the galaxy spectral type 
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obtained from a Principal Component Analysis of the con- 
tinuum and the absorption features. The correlation is far 
better with galaxy spectral types than with morphological 
types. We also found that the extinction derived from the 
emission line B aimer decrement correlates with the galaxy 
spectral type, being larger for early type galaxies. 

But the sample was small: there were only 15 nor- 
mal spirals in the Kennicutt (1992) atlas. In the mean- 
time, Jansen et al. (2000a and 2000b) published a com- 
prehensive photometric and spectrophotometry atlas of 
200 nearby field galaxies (the Nearby Field Galaxy Sur- 
vey, hereinafter NFGS). The objects were selected so as 
to span the full range in galaxy luminosities and mor- 
phological types. As in Kennicutt (1992), the integrated 
galaxy spectra were obtained by drifting a long slit across 
the galaxy. This considerably richer and more appropri- 
ate data base motivated us to revisit the study of SS99, 
in order to see whether the formerly found trends are con- 
firmed, and to refine their empirical interpretation. 

In this paper, we essentially investigate the overall ef- 
fective extinction and metallicity of spiral galaxies through 
an empirical analysis of the properties of their global emis- 
sion lines. We also perform a Principal Component Anal- 
ysis in order to understand the relation amongst the vari- 
ables describing the galaxies. 

As in SS99, our paper deals only with normal galax- 
ies. Therefore, we have rejected from the NFGS sample all 
the galaxies with Seyfert nuclei and nuclear starbursts, as 
well as the galaxies with indication of close interaction or 
strong peculiarity, as indicated in the appendix of Jansen 
et al. (2000a). In the following, the expression NFGS sam- 
ple will refer to the 158 galaxies that have no evident sign 
of peculiarity. 

All the galaxies in our sample are found at redshifts 
lower than 0.05. Their blue absolute magnitudes range be- 
tween -22 and -13. For most of the diagrams investigated 
in this paper, relevant data are available for 100-120 ob- 
jects. 

2. Dust extinction 

The dust content of galaxies is a controversial subject. 
While the first comprehensive analysis of this problem 
(Holmberg 1958, 1975; de Vaucouleurs 1959) have pro- 
posed that spiral disks are largely transparent, more re- 
cent works have suggested that these disks are opaque or 
at least far more obscured than previously thought (cf. 
Valentijn 1990; Disney, Davies & Phillipps, 1989; James 
& Puxlcy 1993). So far, there is still no consensus on 
whether galaxies are optically thick or thin (see Wang & 
Heckman 1996 or Meurer & Seibert 2001 and references 
therein). Actually, the galaxy extinction depends not only 
on the amount of dust and its composition, but also on 
the distribution of dust relative to the light sources. This 
point too has been disputed. For instance, some authors 
suggest a foreground screen dust geometry (Calzetti, Kin- 



ney & Storchi-Bcrgmann 1994, 1996; but see Witt, Thron- 
son & Capuano 1992), while others propose hybrid mod- 
els, with the dust partially distributed in a foreground 
screen and partially concentrated in the star forming re- 
gions (e.g. Chariot & Fall 2000). Indeed, Chariot & Fall 
(2000; see also Chariot & Longhctti 2001) present a model 
for the effects of dust on integrated galaxy spectra that ex- 
plains several observational properties of starburst galax- 
ies, showing that there is a sequence in the overal dust 
content of the galaxies. Radial gradients should also play a 
role, and Giovanelli et al. (1994) and Peletier et al. (1995) 
propose that the dust distribution presents a strong ra- 
dial dependence, with considerable obscuration near the 
centers of galaxies and little at the outer edges (see also 
Nelson, Zaritsky and Cutri 1998). 

In nearby starburst galaxies, there is evidence from 
bidimensional spectroscopy that stars, gas and dust are 
decoupled (Mafz-Apellaniz et al. 1998) and the attenua- 
tion inferred from the Ha /H/3 ratio is typically higher than 
that inferred from the spectral continuum (e.g., Fanelli, 
O'Conncll & Thuan 1988; Calzetti, Kinney & Storchi- 
Bergmann 1994; Mayya & Prahbu 1996; Calzetti 1997). 
Recently, Poggianti et al. (1999) and Poggianti & Wu 
(2000) have shown that the so-called e(a) spectrum can 
be reproduced assuming that in a current starburst the 
dust extinction affects the young stellar populations much 
more than the older stars. This effect has been called se- 
lective dust extinction by Poggianti and collaborators, and 
is a consequence of the fact that a large fraction - but not 
all - of the dust in galaxies is associated with star forma- 
tion regions, absorbing a significant fraction of the light 
emitted by the young stars. We now present an empirical 
evidence that selective extinction is indeed affecting the 
integrated spectra of spiral galaxies. 

2.1. Selective Extinction 

Fig. 1 shows the emission line equivalent width of H/3, 
EW(H/3), as a function of that of Ha, EW(Ha), for the 110 
galaxies from the NFGS sample with available data. The 
shape of the symbols correspond to different morphologi- 
cal classes, as indicated in the figure caption. As explained 
by Jansen et al. (2000b), these equivalent widths were 
extracted from the spectra by placing the limits of the 
measurement window well inside the absorption trough, 
thus minimizing the effect of the underlying stellar ab- 
sorption, which becomes appreciable at small equivalent 
widths. Following the prescription of Jansen, Franx & Fab- 
ricant (2001), we have applied an additional correction for 
residual absorption (1.5 A for Ha and 1 A for H/3). 

The correlation between EW(H/3) and EW(Ha) is im- 
pressively strong and is well represented by a relation of 
the type 

EW{R/3) = A x £W(Ha) (1) 
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Fig. 1. EW(H/3) as a function of EW(Ha) for the normal 
spiral galaxies from the NFGS sample. Different symbols 
correspond to galaxies of different morphological types T. 
Circles: — 5 < T < 0; squares: 1 < T < 3; triangles: 
4 < T < 6; Circles: 7 < T < 10. The total number of the 
objects appearing in the diagram is given in parenthesis. 
Such a presentation is adopted for all the observational 
diagrams in this paper. The EW(H/3) = 0.245 EW(Ha) 
line is shown. 



where A is a constant. Using the ordinary least squares 
bisector estimation for the fitting (Isobe et al. 1990), we 
find A = 0.245±0.007. This correlation is even better than 
the one found by SS99 for the Kennicutt sample. The un- 
certainties in the emission line equivalent widths may be 
appreciable for EW(Hq) or EW(H/3) smaller than, say, 
10 A but the value of the slope A is essentially deter- 
mined by larger equivalent widths and is not sensitive to 
these uncertainties. We have performed a non-parametric 
Kendall's r statistical test (Press et al. 1992) to evalu- 
ate quantitatively the quality of the fitting, for the whole 
sample and for the sub-sample with EW(Ha) > 5A and 
EW(H/3) > lOA. The results are presented in Table 1 and 
reveal that the correlation between EW(Ha) and EW(H/3) 
is indeed highly significant for both samples. For an easy 
consultation, the results of this test for all the observa- 
tional diagrams investigated in this study have been sum- 
marized in Table 1. In this table we present, for each plot 
and sample, the number of data points, the value of the 
non parametric correlation coefficient r, and its two-sided 
significance level, prob. Note that small values for prob 
indicates that the correlation or anti-correlation between 
the variables is significant. 

We found that the residuals in the EW(H/3) vs 
EW(Ha) relation are independent of the galaxy ellipticity 
(therefore inclination) as well as of the galaxy absolute 
magnitude (and any combination of ellipticity and magni- 
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Fig. 2. Ci as a function of the galaxy synthetic colour 
(B — R) derived from the spectra. Symbols have the same 
meaning as in Fig. I. Open symbols correspond to galaxies 
with EW(Ha)< 10 A or EW(H/3)< 5 A. 



tude). They are also independent of galaxy morphological 
type. 

The correlation shown in Fig. 1 has important impli- 
cations concerning the nature of the internal reddening of 
the galaxies, in particular in terms of a link between the 
effective extinction in the continuum and in the emission 
lines. Let us write an emission line equivalent width as a 
ratio between the emission line intensity / and the flux 
F c at the adjacent continuum in the observed spectrum: 
EW — I/F c . The relation between the observed and in- 
trinsic (i.e., extinction corrected) intensities of Ha and H/3 
is 



/(Ha) 



5 X l0-Ci(f(H a )-f(Kf3)) 



(2) 



where S is the intrinsic intensity ratio of these two lines, 
Ci is a measure of the effective extinction at H/3 in the 
region producing the lines, and /(A) is a function that 
gives the wavelength dependence of the extinction. Let C c 
be a measure of the effective amount of extinction at the 
wavelength of H/3 in the region producing the continuum, 
so that the observed and intrinsic fluxes (denoted by F c 
and respectively) at a given wavelength are related by 



F c {\) = F C °(A) x 10- C ^ A ) 



(3) 



Eq. 1 implies that the effective extinction at H/3 is related 
to the ratio of the observed continuum fluxes at Ha and 
H/3 by 



Ci 



-1 



/(Ha) - /(H/3) 



log 



1 F c (Ha) \ 
A5F C (K0)J 



(4) 
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and that the difference between the line and continuum 
effective extinctions is related to the ratio of the intrinsic 
continuum fluxes by 



Ci - C c 



■1 



/(Ha)-/(H/3) g \A5F?(m) 



(5) 



Since F® (Ha) / F® (H/3) depends only on the stellar pop- 
ulations of the galaxy, Eq. 5 is indeed an empirical pre- 
scription for the relation between selective dust extinction 
and the galaxy stellar populations. 

Although we have no way to determine the value of 
C c from the available data, some implications of Eq. 5 
may be inferred. First, the difference between C'i and 
C c is essentially dependent on the intrinsic colour of the 
galaxy and is larger for intrinsically redder galaxies. Sec- 
ond, if Ci > C c , a lower limit for the intrinsic con- 
tinuum flux ratio at the wavelengths of Ha and H/3 is 
F c °(Ha)/F c °(H/3) = AS. Assuming S = 2.9 (the case B 
recombination value corresponding to a temperature of 
9000K), and using the Seaton (1979) extinction law, we 
find that F c °(Ha)/F c °(H/?) in spiral galaxies should be 
larger than 0.70. 

We have used the software PEGASE.2 (Fioc & Rocca- 
Volmerange 1997) with a Salpeter IMF and default val- 
ues to compute this ratio for a few simple model galaxy 
spectra. For an instantaneous burst 15 Gyr old we ob- 
tain F c °(Ha)/F c °(H/3) = 0.84 and Eq. 5 implies d - C c = 
0.21, whereas for a model with constant star formation 
rate and same age this ratio is 0.75 and Eq. 5 leads to 
Ci — C c — 0.07. Additional clues on C c may be obtained 
from spectral synthesis methods that include explicitly in 
the analysis the effects of extinction on the spectra (e.g., 
Cid Fernandes et al. 2001). 

2. 2. The extinction of the nebular light 

We obtained C; from the observed Ha/H/3 ratio for the 
110 normal galaxies of the NFGS sample with available 
data. 

Fig. 2 plots Ci as a function of the galaxy color (B—R). 
The open symbols correspond to galaxies with EW(Ha)< 
10 A or EW(H/3)< 5 A for which the value of C t tends to 
be more uncertain. A very good correlation is seen, with 
d larger for redder objects (the correlation coefficient t 
shown in Table 1 is larger than 0.6). This behavior is in- 
deed expected from Eq. 4. As already shown by Jansen, 
Franx & Fabricant (2000), there is a correlation between 
d and the galaxy blue absolute magnitude, although the 
dispersion is much larger than for the color. C; also cor- 
relates with the morphological types, decreasing towards 
later types, confirming the trend found in SS99. The max- 
imum extinction is observed for galaxies Sb, but the vari- 
ance is large. On the other hand, there is no correlation 
whatsoever with the galaxy apparent ellipticity, implying 
that at least part of the obscuring dust is concentrated 
close to the emitting sources. 
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Fig. 3. Ci as a function of [O in] A5007,4959 + [O n] 
A3726,3729) / H/3 corrected for reddening. Same symbols 
as in Fig. 2. 

We now examine whether there is a relation between 
Ci and the overall metallicity of the galaxies. So far, there 
have been contradictory claims in this respect. Zaritsky, 
Kennicutt & Huchra (1994) found no evidence for a sys- 
tematic dependence between reddening and abundance in 
a sample of 39 disk galaxies. In their study, the character- 
istic metallicity of the galaxies was derived from the value 
of the metallicity indicator ([O ill] A5007,4959 + [O n] 
A3726,3729) / H/3 (hereinafter 23 ) (Pagel et al. 1979, 
McGaugh 1991) at the isophotal radius, interpolated from 
the values observed in giant H II regions. In other contexts, 
however, it has been suggested that the extinction derived 
from the Balmer decrement might be related to metallic- 
ity. For example, Campbell, Terlevich & Melnick (1986) 
found a loose correlation between the extinction and the 
oxygen abundance in a sample of H n galaxies. Also, in- 
vidivual giant H II regions in spiral galaxies exhibit a de- 
crease of extinction with galactocentric radius (Kennicutt 
& Garnett 1996; Van Zee et al. 1998a), but this is not a 
general rule (Martin & Roy, 1992 and references therein). 

Fig. 3 shows Ci as a function of the usual metallicity 
indicator 023- The meaning of the symbols is the same 
as in Fig. 2. When considering only the filled symbols, 
there is a weak correlation in the sense that 23 tends 
to be smaller for larger values of Cj. The correlation vir- 
tually disappears when including the open symbols, for 
which the extinction and reddening corrected line ratios 
tend to be more uncertain. These results are confirmed by 
Kendall's r statistical test; see Table 1. We know that the 
relation between O23 and O/H is double valued, with a 
maximum in <3 23 occuring at log O/H + 12 of about 8.2. 
This is illustrated in Fig. 4, which shows O23 versus O/H 
as computed for a sequence of simple constant density 
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Fig. 4. [O in] A5007,4959 + [O n] A3726,3729) / H/3 as 
a function of O/H for a series of photoionization models 
representative of giant H n regions (total number of ioniz- 
ing photons 3 10 51 ph s _1 , effective temperature 45 000 K, 
density 100 cm~ 3 and filling factor 0.03). 



models roughly representative of giant H n regions, using 
the photoionization code PHOTO described in Stasihska 
& Lcithcrcr (1996). Because radial abundance gradients 
are known to exist in spiral galaxies from the study of 



their bright giant H II regions (Vilas-Costas & Edmunds, 
1992, Zaritsky, Kennicutt & Huchra 1994), it is conceiv- 
able that, depending on the abundance gradients and the 
relative weights of the regions of various abundances in the 
integrated emission line fluxes, O23 might not be a a good 
indicator of the galaxy overall metallicity. It is important 
to note that metallicity here refers only to the star for- 
mation regions. Kobulnicky, Kennicutt & Pizagno (1999) 
have addressed this problem by simulating the global O23 
of 22 spiral galaxies, integrating the observed values for 
several bright individual H 11 regions over the Ha radial 
profiles from CCD images. They found that the oxygen 
abundance derived from the simulated integrated O23 was 
identical, within errors, to the characteristic abundance of 
galaxies defined by Zaritsky, Kennicutt & Huchra (1994) 
as the oxygen abundance computed at a certain character- 
istic galactocentric radius. They concluded that the beam 
smearing effect from sampling large number of galaxies, 
even in the presence of strong abundance gradients, has 
a small effect on characterizing the mean abundances of 
galaxies. They mention however, that it would be useful to 
test this conclusion directly using actual integrated spec- 
tra of galaxies. Unfortunately, none of the galaxies in their 
sample belongs to the NFGS sample. Our conclusion, us- 
ing the NFGS sample is purely empirical, and is that O23 
docs not seem to be a useful indicator of metallicity in 
integrated spectra. 

Another indirect empirical metallicity indicator has 
been proposed more recently (van Zee et al. 1998a), which 
has the double advantage of not being plagued by the dou- 
ble value problem and of being independent of reddening. 
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Fig. 7. Ci as a function of [N n] A6584/[0 n] A3727 cor- 
rected for reddening. Same symbols as in Fig. 2. 
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Fig. 6. Ci as a function of [N n] A6584/Ha corrected for 
reddening. Same symbols as in Fig. 2. 



This is the [N n] A6584/Ha ratio. Fig. 5, obtained from the 
same models as Fig. 4 (these models assume that N/O a 
(O/H) - 5 ) shows that [N n] A6584/Ha increases with O/H 
up log O/H + 12 around 8.8. There are two reasons for 
such a behaviour. One is that [N n] A6584/Ha is less de- 
pendent on the electron temperature than O23, so that it 
is not so much affected by the decrease in electron tem- 
perature provoked by an increase in metallicity. The other 
reason is that N/H increases with O/H. The study of giant 
H II regions in spiral galaxies (van Zee, Salzer & Haynes 
1998b) indicates that at log O/H + 12 > 8.5, N/O is pro- 



portional to O/H and not to (O/H) 5 , so that the effect 
is even stronger than in the models shown in Fig. 5. 

Fig. 6 shows Ci as a function of [N 11] A6584/Ha for 
the normal galaxies in the NFGS sample. There is a very 
clear correlation! Kendall's r is 0.6 and the small value 
of prob indicates that the correlation is highly significant. 
As a matter of fact, it is even surprising that such a good 
correlation should exist. Indeed, the [N 11] A6584/Ha ra- 
tio obviously strongly depends on the ionization param- 
eter of the H II regions. While giant H II regions seem 
to follow a very narrow sequence in ionization param- 
eters (McCall, Rybski & Shields 1985; Dopita & Evans 
1986; Stasinska, Schaerer & Leitherer 2001), one expects 
in integrated spectra of spiral galaxies a non negligible 
contribution from a diffuse medium, possibly ionized by 
photons leaking out from bona fide giant H II regions. 
Indeed, about 20 to 50% of the H-Balmer line emission 
in nearby spiral galaxies arises from this diffuse medium 
(e.g., Greenawalt, Walterbos & Braun 1997; Collins et al. 
2000; Rossa & Dettmar 2000; Zurita, Rozas & Beckman 
2000; Rand 2000). The [N 11] A6584/Ha ratios observed 
in this diffuse medium are typically a factor 2 larger than 
in the nearby giant H 11 regions, which implies that in an 
integrated spectrum, the effect of the diffuse medium at 
a given galactocentric distance (and therefore metallicity) 
is to increase the [N 11] A6584/Ha ratio by about 30% at 
most. Such an increment is relatively modest compared to 
the metallicity dependence of this ratio and explains why 
[N 11] A6584/Ha appears to be such a good indicator of 
the overall metallicity of spiral galaxies. We should em- 
phasize that, nevertheless, this indicator cannot be easily 
translated into a value of O/H, in the presence of abun- 
dance gradients and integration over the complex inter- 
stellar medium of spiral galaxies. For example, the weight 
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of the regions of mctallicities around solar is expected to 
be important, as seen from Fig. 5. One should certainly 
refrain from using the calibration that has been proposed 
by van Zee at al. (1998a), which is valid only in the context 
of giant H n regions. 

Another potential empirical indicator of metallicity of 
the emission line regions which is independent of the ion- 
ization parameter, does not suffer from the double value 
problem and is less affected by metallicity gradients is 
[N n] A6584/[0 n] A3727 (Dopita et al. 2000). It is how- 
ever strongly dependent on the reddening correction. Fig. 
7 shows d as a function of [N II] A6584/[0 II] A3727 (the 
open symbols have the same meaning as in Figs. 2, 3 and 
6), and here again we clearly see a correlation. It is more 
dispersed than the C; versus [N n] A6584/Ha relation, pos- 
sibly because of errors in the reddening correction. Indeed, 
the correlation is much better when considering only the 
filled symbols. This is confirmed by the results of Kendall's 
t test presented in Table 1. 

From the above considerations, our conclusion is that, 
indeed, Ci correlates with the global metallicity of the 
emission line regions in spiral galaxies, at variance with 
the conclusions of Zaritsky, Kennicutt & Huchra (1994). 
It should be stressed that C; does not necessarily measure 
the true extinction of the total light emitted in the Balmer 
lines. In the case of a patchy dust distribution, it measures, 
so to say, the effective extinction of those regions that are 
not too strongly obscured, so that the interpretation of 
our result is by no means straightforward. 



3. The emission line sequence 

In SS99, we showed that, when ordering the normal galax- 
ies from the Kennicutt sample by galaxy spectral type, 
emission line equivalent widths and emission line ratios 
from the integrated galaxy spectra formed a very well 
defined sequence (not necessarily monotonic!). We inter- 
preted this as an evidence for a close link between the 
high mass end of the stellar populations, which drive the 
emission lines, and the population of the stars that are 
responsible for the optical continuum, which determines 
the galaxy spectral type. The scatter in the properties of 
individual H II regions observed in spiral galaxies (Kenni- 
cutt & Garnctt 1996, Roy & Walsh 1997) is likely to be 
smoothed out in integrated galaxy spectra if the galaxies 
are sufficiently massive to contain a large number of H II 
regions. The sequence found by SS99 was so impressive, 
that it calls for further investigation and the NFGS sample 
is adequately large for such a purpose. Unfortunately, we 
cannot determine the galaxy spectral type from the pub- 
lished data, since the spectra are not yet publicly available 
in digital form. As a proxy of the galaxy spectral type, 
we can use galaxy colours, with the advantage that they 
are much easier to obtain than spectral types. We found 
empirically that the best colour to use was the synthetic 



(B — R) derived from the galaxy spectra as explained by 
Janscn ct al. (2000b). 

3.1. Equivalent widths and star formation 

In Fig. 8, we show the equivalent widths of Ha, [O n] 
A3727, [O in] A5007 and [N n] A6584 as a function of 
(B — R). The shape of the symbols represent the galaxy 
morphological types, as in the previous observational dia- 
grams. Clearly, the equivalent widths of Ha, [O n] A3727 
and [O in] A5007 are strongly correlated with (B — R), be- 
ing larger for bluer galaxies. However, there seems to be 
more dispersion than what was suggested by the study of 
SS99. One possible reason is that the color (B — R) is not 
as good a parameter as the spectral type for ranking galax- 
ies. For comparison, we show in Fig. 9 the same equivalent 
widths but as a function of the galaxy morphological type, 
T. The correlations of equivalent widths with T are seen 
to be much weaker than with (B — R). These conclusions 
are confirmed by the statistical test summarized in Table 
1. 

Emission lines, principally Ha and [O n] A3727, are 
indicators of current star formation (Barbaro & Poggianti 
1997, Kennicutt 1998, Jansen, Franx & Fabricant 2001). 
The calibration in terms of star formation rates how- 
ever strongly depends on several factors: star formation 
law, metallicity, selective dust extinction, etc. Chariot & 
Longhctti (2001) argue for the necessity of using the [O n] 
A3727, [O in] A5007 Ha and H/3 lines at the same time to 
better constrain the problem. Qualitatively, though, Fig. 8 
(together with Fig. 1) shows that any of these parameters 
is likely to be a first order ranking indicator of integrated 
star formation rates in spiral galaxies. This is of course 
not the case of [N n] A6584, whose equivalent width be- 
haves in a completely different manner. Indeed, the strong 
drop in nitrogen abundance for the bluest (less chemically 
evolved) galaxies overtakes the effects of increasing present 
star formation rate and increasing electron temperature. 

As shown by Jansen et al. (2000b), equivalent widths 
are also loosely correlated with galaxy absolute magni- 
tudes, tending to be larger for more luminous galaxies. 
However, at a given galaxy colour, we found no com- 
pelling evidence of brighter galaxies having larger equiva- 
lent widths. 

3.2. Line ratios and metallicity 

Fig. 10 shows various emission line ratios as a function 
of (B — R). The line ratios have been dereddened using 
the value of C; computed above. Since the reddening cor- 
rection becomes more uncertain for galaxies with small 
EW(Ha) and EW(H/3), as in Figs. 2, 3, 6 & 7 we marked 
by filled symbols those objects that have EW(Ha)> 10 A 
and EW(H/3)> 5 A and thus the most reliable reddening 
corrections. Clear trends are seen for [O in] A5007/[O n] 
A3727, [N n] A6584/[0 n] A3727 and [N n] A6584/Ha as 
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a function of (B — R). As for the equivalent widths, we 
find the trends much clearer with (B — R) than with mor- 
phological type. For [O ill] A5007/[O n] A3727 and [N n] 
A6584/[0 n] A3727, which are sensitive to reddening, one 
can see that the open symbols prolongate the trends de- 
fined by the filled symbols, but with much higher disper- 
sion. In the case of O23, there is a weak trend with (B — R) 
when considering only the filled symbols, and practically 
no trend when considering the entire data set. The sta- 
tistical significance of these trends are presented in Table 
1. 



The fact that [O 111] A5007/[O 11] A3727 decreases as 
the galaxy color gets redder is not necessarily only due 
to a lowering of the average ionization parameter or the 
average spectral hardness of the exciting stars (as was sug- 
gested by SS99). Indeed, an increase in overall metallicity, 
by enhancing the temperature drop in the central zones of 
H II regions with metallicities larger than solar, also con- 
tributes to lower the [O ill] A5007/[O 11] A3727 ratio (see 
Stasinska, Schaerer & Leitherer 2001). We find no trend 
of [S 11] /Ha with colour, suggesting that the effect of ion- 



G. Stasiriska & L. Sodre Jr.: The reddening and metallicity sequence in spiral galaxies 



O 
O 



1 1 1 1 1 1 1 1 1 1 

: <-) . . . 


■ i i i i 

(we), : 


■ ii 

■ i : 


: * * * - 

A * * 7 

★ * 




* * 


■ ■ ■ 




■ 

■ ■ 

• a ■ 

- • ■ 
• 




■ 

• ■ 

i i i i i i i i i i i 





5 
T 



10 



o 



o 



n 1 1 1 1 1 1 1 r- 



(b) 



■> — i — ' — r~- 
(123)* i 



* i * 



t A * 

• • * * * 



■ s * * 

: a ! 4 * : * $ 

■ | t A * * * * " 



_i i ? i i t ■ ■ l ~ 



5 
T 



10 



o 
o 



= ° 
o "~ 



I 1 1 1 1 I 1 1 1 1 I 1 

: (c) 


1 1 1 1 - 

( 81) : 


■ 

• 


* 


A 
A 


* * * ★ - 


• * 

a ! . 

• * 


★ 

* * * - 
* ★ - 

★ 


■ i • . 


★ ★ 
★ 


■ • i 


★ 


A A 

— 5 — i — i — i — i — | — i — i — i — i — | — 1_ 


i i i i 



-5 







10 



T 



o 



F — i — 1 — 1 — 1 — 1 — r 



^= ° 

LjJ 



i 1 1 r 



(d) 



(116) 



A i - 

■ * A 

■ A A a 

I " 4 * * * 

A * * * $ 



* * * * 

* A A * , 

* 



★ ★ 



— i i i A i i i i_ 



-5 







10 



T 
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ization parameter along the colour sequence may not be 
overwhelming. 

As mentioned in Sect. 2, the O23 ratio is not the best 
indicator of metallicity in integrated spectra of galaxies 
where abundance gradients are known to occur. Being non 
monotonic with respect to the metallicities of individual 
H II regions, when integrated over an entire galaxy it be- 
comes impossible to interpret without additional informa- 
tion on the radial galactic gradients in abundances and 
star formation rates. In addition, it is strongly dependent 
on the reddening correction, which becomes uncertain in 
metal rich galaxies with weak lines. 



On the other hand, [N 11] A6584/[0 11] A3727 is a good 
indirect indicator of metallicity in extragalactic H 11 re- 
gions (Dopita et al. 2000) and is expected to be a fair in- 
dicator of overall metallicity in integrated galaxy spectra 
(see above). The trend seen in Fig. 10c is produced by an 
increase of N/O as galaxy colour gets redder, indirectly 
indicating an increase in metallicity. As a consequence, 
the conspicuous trend of [N 11] A6584/Ha increasing with 
(B — R) (Fig. lOd), is likely due to an increase in metal- 
licity. As mentioned before, the [N 11] A6584/Ha ratio has 
the advantage of being reddening independent, and it can 
be used with confidence even for galaxies with weak emis- 
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sion lines. It clearly indicates that the sequence between 
metallicity and (B — R) colour is quite narrow even for 
galaxies redder than (B — R)=1.2 (which is not appar- 
ent when using the strongly reddening-dependent [N ii] 
A6584/[0 n] A3727 ratio). 

Numerous studies have claimed the existence of a 
metallicity-luminosity relation in a variety of classes of 
galaxies: dynamically hot galaxies, i.e. ellipticals, bulges, 
and dwarf spheroidals, dwarf H n galaxies, irregular galax- 
ies and spirals (Bender et al. 1993; Vigroux, Stasiriska & 
Comte 1987; Skillman, Kcnnicutt & Hodge 1989; Richer 
& McCall 1995; Zaritsky, Kennicutt & Huchra 1994). In 



Fig. 11, we show the various metallicity indexes, O23, [N Ii] 
A6584/Ha and [N 11] A6584/[0 11] A3727 as a function of 
the total absolute blue magnitude Mb ■ We do find a strong 
correlation between Mb and the metallicity indexes [N 11] 
A6584/Ha and [N 11] A6584/[0 11] A3727. The correlation 
of the O23 index with Mb is statistically significant (see 
Table 1), but at a rather low level. It is interesting to recall 
that Zaritsky, Kennicutt & Huchra (1994), by using this 
index to infer the characteristic metallicity of spiral galax- 
ies, did find a trend with galaxy luminosity. But the index 
actually referred to the observed ratio at a certain galac- 
tic radius, and not to the integrated light from the galaxy. 
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Our study has shown how misleading can be the use of O23 
in integrated spectra of galaxies in the presence of abun- 
dance gradients. On the other hand, [N 11] A6584/Ha and 
[N 11] A6584/[0 11] A3727 show clear trends with galaxy 
absolute magnitude, confirming that indeed, there is a re- 
lation in spiral galaxies between the overall metallicity of 
the star forming regions and the galaxy luminosity. 

4. Principal Component Analysis 

Principal Component Analysis (PCA) is a useful tool for 
the examination of multidimensional data and the iden- 
tification of underlying variables that may be responsible 



for the variance in a data set (Murtagh & Heck 1987, 
Fukunaga 1990). PCA allows to define a new orthonormal 
reference system in a parameter space, with basis- vectors 
(the principal components) spanning directions of maxi- 
mum variance. 

Here we apply PCA to a set of variables describing 
different aspects of the galaxies in the NFGS sample. The 
7 variables selected for this exercise are: the blue absolute 
magnitude Mb] the color (B — R); the numerical value 
of the morphological type T; the galaxy surface bright- 
ness at the effective radius /i B ; EW(Ha), which is related 
to the star formation rate; the metallicity indicator [N 11] 
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A6584/Ha; and Hcv/H/3, related to the emission-line ex- 
tinction. The number of normal galaxies with complete 
information that enters in the analysis is 102. 

Since the input variables are of different nature, it is 
appropriate to apply PCA on the correlation matrix, that 
is, the input variables were transformed: they had their 
mean value subtracted and, after, they were divided by 
their standard deviation. Hence the analysis is made on 
a set of transformed variables with zero mean and unity 
variance. The result indicates that the first principal com- 
ponent explains alone 48% of the sample variance. The 
second and third components are responsible for 19% and 
12% of the sample variance, respectively. 

Although each principal component is a linear combi- 
nation of all the input variables, it is interesting to verify 
whether some variables are well correlated with the first 
components because, in this case, they are responsible for 
a significant fraction of the sample variance. We have then 
computed, for each variable and principal component, the 
non-parametric Kendall's r coefficient of rank correlation 
and the respective probability of absence of correlation 
between these quantities. We found that the variable that 
best correlates with the first component is the metallic- 
ity indicator [N n] A6584/Ha; the variables (B — R), and 
Ha/H/3 also correlate well with the first component; the 
morphological type is correlated with the first principal 
component, but to a lesser extent. On the other side, only 
fj, e B is strongly correlated with the second principal com- 
ponent. The absolute magnitude correlates best with the 
third principal component. Hence, this analysis indicates 
that the variance in the parameter space considered here is 
produced, in first place, by the metallicity and parameters 
related to it or to the stellar populations, and, in second 
place, by the surface brightness, which is related to the dy- 
namical history of the galaxies. The absolute magnitude, 
related to the mass of the galaxy, surprisingly comes only 
in the third place. 

5. Conclusions 

We have explored the emission line trends from the inte- 
grated spectra of normal spiral galaxies, in order to inves- 
tigate the relationships between dust extinction, metallic- 
ity and some macroscopic properties of spiral galaxies. A 
similar aim was followed by the study by Zaritsky, Ken- 
nicutt & Huchra (1994) on a sample of 39 disk galax- 
ies, which however was not based on integrated spectra of 
galaxies, but on the study of several giant H II regions in 
each galaxy. Our approach, using integrated spectra, is a 
priori likely to give more robust answers. It also has the 
advantage of being readily applicable to studies of galaxies 
at high redshift. We benefited from data from the Nearby 
Field Galaxy Survey (Jansen et al. 2000a and b), which 
considerably increased the size of the so far available sam- 
ples for systematic studies of the global properties of spiral 
galaxies. 



Our main results are the following. 

We found that the correlation between EW(Ha) and 
EW(H/3) is extremely strong. This implies that Ci -C c , 
the difference between the effective extinction of the stellar 
and the nebular light, depends only on the intrinsic colours 
of the galaxies, being larger for redder galaxies. 

The nebular extinction as derived from the Balmcr 
decrement is found to correlate essentially with the effec- 
tive metallicity of the emission line regions. This finding 
required using an appropriate metallicity indicator. The 
usual O23 index does not work for integrated spectra of 
galaxies, due to strong metallicity gradients. On the other 
hand, [N 11] A6584/Ha and [N 11] A6584/[0 11] A3727 seem 
to be adequate empirical metallicity indicators, the first 
one having the advantage of being reddening independent 
and being applicable also for galaxies with weak emission 
lines. The calibration of these indicators would require 
a modelling of the integrated spectra of galaxies, in the 
vein of the works of Chariot & Longhetti (2001) or Moy, 
Rocca-Volmerange & Fioc (2001) but including the effects 
of radial gradients in the metallicity and star formation 
rates. 

We find that the overall metallicity of the emission line 
regions, as inferred from adequate indicators, is strongly 
correlated with galaxy colours. The correlation is much 
stronger than with morphological types. This suggests 
that in normal spiral galaxies, the global metallicity of 
the star forming regions and the old stellar populations 
are closely linked together, while morphological types are 
influenced by additional factors in the galaxy histories. 

In order to better evaluate the relation between the dif- 
ferent parameters describing the galaxies, we performed a 
PCA analysis on a 7-D parameter space. We found that 
the variance in the considered parameter space is pro- 
duced, in first place, by the metallicity and parameters 
linked to the stellar populations, and, in second place, by 
the surface brightness, which is linked to the dynamical 
history of the galaxies. The absolute magnitude, related 
to the mass of the galaxy, comes only in the third place. 

The strong empirical evidences provided by this ex- 
ploratory work bring important constraints that will have 
to be reproduced by models of spiral galaxy formation and 
evolution, in order to understand what drives the spiral 
galaxy sequence. 
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